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Abstract 

We study the effect of massive isocurvaton on density perturbations in quasi-single field 
inflation models, when the mass of the isocurvaton M becomes larger than the order of 
the Hubble parameter H . We analytically compute the correction to the power spectrum, 
leading order in coupling but exact for all values of mass. This verifies the previous numerical 
results for the range < M < 3H/2 and shows that, in the large mass limit, the correction 
is of order H 2 /M 2 . 



1 Introduction 



Massive fields are present during inflation, as a consequence of realistic model-building and 
UV completion. In terms of the low energy effective field theory, the inflaton travels along a 
trajectory determined by the minima of some complicated potential landscape spanned by a 
large number of fields. The scalar fields orthogonal to this trajectory, the isocurvatons, have 
a rich spectrum. In inflation scenarios from string theory and supergravity, the lightest of 
them start from order H, the Hubble parameter during inflation. This value is determined 
by the coupling of the scalar fields to the spacetime curvature and stabilized by supersym- 
metry. Lighter values require fine-tuning or extra symmetry. Those inflation models with 
isocurvaton mass of order 0(H) are called the quasi-single field inflation models [l}|5]. The 
simplest version of such models has been used to show the distinctive effect of these massive 
isocurvatons on the primordial non-Gaussianities. In this paper, we would like to consider 
the model for a different purpose. We will investigate how exactly the massive field decouples 
when its mass becomes larger than 0(H). 

In addition to the isocurvatons with mass of order 0(H), there are certainly more fields 
that have masses larger than 0(H). Intuitively the quantum fluctuations of such fields should 
be less important to the low energy effective theory description. The purpose of this paper 
is to make this statement more precise and quantitative in the context of the quasi-single 
field inflation, which represents a quite general class of models in model-building. The effects 
of heavy fields have been studied in different models previously [6j[T] . It has been argued 
that their corrections to the power spectrum are proportional to powers of H/M for large 
M/H, and the specific power is model-dependent. In the simplest model of quasi-single field 
inflation studied in [lj[2], the massive field couples when the inflaton trajectory turns, and 
the coupling appears in the kinetic term. We would like to compute the effect of this massive 
field on power spectrum for all values of M. Generally speaking, obtaining such a correction 
in an analytically closed form from the first principles is difficult, and such an example was 
not available so far. This is because the in-in formalism involves multiple integrations of 
the product of the mode functions, and the mode function for the massive field contains a 
special function. It turns out that the model example studied in [TJ|2) is just simple enough 
for us to achieve this goal here, thus providing a solid, and reasonably general case on how 
precisely we mean that massive fields can be integrated out. 

Before the detailed computation, we may have two different guesses what the final result 
should be qualitatively. The first possibility is by analogy to the thermal field theory, in which 
the contributions of massive states to correlation functions are exponentially suppressed by a 
Boltzmann factor if the mass is much higher than the temperature. Here in de Sitter space, 
we have a Gibbons-Hawking temperature Tqh — H/2n, so the corresponding Boltzmann 
factor is e ~ 27TM / H _ The second is by analogy to the scattering process in particle physics, in 
which the contribution of massive states to the scattering amplitude is power-law suppressed 
if the mass is much higher than the energy scale of the process. Interestingly as we will see, 
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both types of contributions appear in different terms in the in-in formalism, and the latter 
is more general. Consequently, for large M/H, the power-law suppression dominates. 



2 The model and a summary of the results 

For convenience we review the model and summarize the main results in this section. 

We consider the constant turn case of a two-field quasi-single field inflation model [Tlf2] , 
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d^xsj—g 



-\(R + crfg^d^e - \g^d^ad v a - V ST {0) ~ V(a) 



(2.1) 



This describes an inflaton that is moving around a circle in field space with radius R, and at 
the same time rolling down the potential. The 9 is the angular coordinate, and a is the radial 
coordinate. In the angular direction, the potential is a slow-roll potential V ST (9); while in the 
radial direction the potential is lifted around the effective minimum <7o so the isocurvaton is 
massive. The massive field not only determines the bending of the trajectory classically, but 
also influences the fluctuations of the inflaton field quantum-mechanically. In the infinite 
mass limit M — > oo, this Lagrangian reduces to a single field model, with radius fixed at 
R = R + o"o and the inflaton field being R9. The purpose of this paper is to demonstrate 
how this limit is reached and what the dependence of power spectrum on the mass M is for 
a finite M. 

After perturbing the two fields, we obtain the following kinematic Hamiltonian density 
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and the interaction Hamiltonian density 



(2.2) 



(2.3) 



0: 



We have defined M to be the effective mass appearing in the Hamiltonian, M 2 = l /// (o"o)+3^ ( 
and 9q is the constant turning angular velocity. We are only interested in the quadratic part 
of the interaction Hamiltonian in this paper. The label "J" denotes that the fields are in the 
interaction picture. 
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Figure 1: The correction to the leading power spectrum. 



The quadratic interaction Hamiltonian introduces a coupling between the massive isocur- 
vaton and the inflaton, and the strength of the coupling is constant and proportional to 6*0. 
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Due to the constant turn assumption, we are excluding various phenomena associated with 
sharp features in turning trajectory (8 13 . The leading order correction to the power spec- 
trum comes from the interaction illustrated by the diagram Fig. [TJ and the corresponding 
terms in the in-in formalism (see (14)|15) for review) are 



t rt 

dh I dti(0|#/(ti) 56j iZj(ii)|0) 

to J to 

2 Rc [ dh I' dt 2 {U\89 2 I H I {t 1 )H I {t 2 )\0) 

to J to 



(2.4) 



We assume the Bunch-Davies vacuum for both mode functions, 58j and 5ai. The final power 
spectrum is given by p] 
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1 + 8C 
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H 



(2.5) 



In [2], the coefficient C is computed numerically for M/H < 3/2. In this work we will 
compute this coefficient for all M analytically. The coefficient C contains two contributions 
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(2.6) 
(2.7) 

l - , (2.8) 



conformal time. is the Hankel function of the first kind. The parameter \x is defined to 
be 



coming from the two terms in (2.4) respectively. The variable x^ = —ik\Ti where r» is the 



/x = 



M 2 
iJ 2 " 
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(2.9) 



It can be either real or imaginary. For real /i, the isocurvaton is heavier, M > 3H/2, which 
is the main interest of this paper. For imaginary fi, this is related to the parameter v in (Tl|2) 
by \x = —iv, and the isocurvaton has a mass < M < 3H/2. Also note that, for the latter 



case, the factors e ^ in (2.7) and (2.8) disappear 



The first term is a product of two integrals without time-ordering, one from the inter- 
acting vacuum and another from its complex conjugate. It turns out that in the large /i 
limit, this term contains a Boltzmann factor C\ — > n 2 e^ 2nM ^ H . The second term involves 
a time-ordering double integral, similar to the form of scattering amplitude in the in-out 
formalism. In the large /x limit, the suppression turns out to be power-law C 2 — > H 2 /(AM 2 ). 
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We expect the power-law to be the more generic behavior in the in-in formalism. Therefore 
overall C — > H 2 / (AM 2 ). This verifies the result from the semi-classical approximation in the 
large M limit (l6||8). 

The complete results with arbitrary fi (i.e. arbitrary isocurvaton mass) are 
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4 COsh (7T/i) 



(2.10) 
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where ip^ is the polygamma function 
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(2.12) 



We plot the coefficient C in Fig. [2} 
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Figure 2: The final result C (solid lines), C\ and Ci (dashed lines) as a function of M/H for 
both M ~ H and M ^> H. The dots are numerical results. 

If we are only interested in the behavior in the large M/H limit, a simpler IR expansion 
method may be used. In this method, we replace the UV part of the mode function of the 
massive field by its IR approximation, which is actually valid between the super-Hubble 
scale kr — > and the sub-Hubble scale kr ~ —^JM/H. This produces the correct answer 
if the difference between the replaced part and the original one is much smaller than the 
final result. We will also show that the dominant contribution from these heavy fields comes 
from between of order ln/z e-folds before the Hubble crossing point kr = —1, and a few 
e-folds after Hubble crossing. Therefore, the main results of this paper also apply to the 
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non-const ant-t urn long as the variation of the angular velocity and mass is not very 

fast, 6/ {HQ) < 1 and M/{HM) < 1. 

The rest of the paper is organized as follows: In the following two sections, Sec. [3] and [4], 
we give the computational details and tricks for the two integrals. In Sec. |5j we analytically 
continue the final result to imaginary fi, and show that they match the expected properties 
and the numerical results in [2] . In Sec. [6j we give a simpler IR expansion method that can 
be used to reproduce the analytical result in the large M/H limit. In Sec. [7j we investigate 
where the dominant contribution of the massive field comes from during its evolution by 
examining the integrals as a function of an IR cutoff x c . In the Appendix, we describe a 
numerical method for checking the integrals. 



3 The integral without time-ordering 

We start with the simpler integral, 
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(3.1) 



The indefinite integration inside the absolute value takes the form 



dx! x~ 1/2 H^{ Xl )e tXl 
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— i^l rg2 

* ^ r (-^) 2F2 + ifi, I + ifi; ^ + ifi, 1 + 2z>; 2zxi j . 



(3.2) 



where 2^*2(01, ^2! 61, 6 2 ; -2) is the hyper-geometric function defined as 

ai{ai + 1) • • • {ax +n — 1) a 2 (a 2 + 1) • • • (a 2 + n — 1) z n 
61(61 + 1) ■ ■ ■ (61 + n - 1) 6 2 (6 2 + 1) ■ ■ ■ (6 2 + n - 1) n! 



2F 2 {a 1 ,a2]b 1 ,b 2 ;z) = ^ 



(3.3) 



n=0 



The Xi = limit of equation (3.2) vanishes. In the x\ — > +00 limit, using the asymptotic 
behavior 

2 F 2 {a 7 a;b l ,b 2 ;z) ^ T MM e ^^ 



r(6i)r(6 2 ) 



T(a)r(6i — a)T(6 2 — a 



T{a) 

i-z)~ a (ln(-z) - ^(61 - a) - ^(6 2 - a) - ifi(a) - 2 7 ) (3.4) 
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where 



d\nT{z) 

dz 



(3.5) 



Eq. (3.2) takes the form 



— i)e * 



(3.6) 



COsh(7T/i) 

where the oscillatory term e 2lXl at infinity X\ — > oo is take to be zero by an ie-description. 



Taking the absolute value and plugging it in (|3.1|), one finds 

Ci 
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4 COsh (7T/i) 



(3.7) 



A plot of Ci as a function of real \x is shown in Fig. [3} We also numerically check the 
integration and plot it in the same figure. The method of numerical calculation is described 
in Appendix [A| 

In the large M/H limit, the above result goes as 



lim d = vrV 2 ^ . 

fl— >co 



(3.8) 



Interestingly, the suppression in C\ is a Boltzmann factor. However, as we shall show in 
the next section, this Boltzmann factor is not the dominate contribution to the two-point 
function. 



4 The integral with time-ordering 

The integral with the time-ordering and double integration is the more difficult one, 
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Figure 3: Comparison of analytical and numerical results. The left, middle, right panels 
are analytical and numerical plots for C%, C2 and C = C% + C2 respectively. It is clear from 
the plot that C\ is decaying exponentially as a function of /i. However the decay of C2 is not 
as fast as exponential. 



The inner layer of (|4.1|) can be integrated out as 
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(4.3) 



The first term above is the x — > 00 limit of the indefinite integral, and the other two terms 
come from the x — > x\ side. 

The difficult part is the integration in the second layer, 



/■oo 

/ dx x x- 1/2 H^{x 1 )e- ix ^X{x 1 ) . 
Jo 



(4.4) 



With both the special function 2F2 and in the integrand, the indefinite integral cannot 
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be done by Mathematica 8. We use a trick of resummation. We expand the part of the 
integrand in the second layer, x\ x ^ 2 (xi)e~ tXl , in a series expansion in bases of x n± ' 1 ^ 1 /' 2 . 
With the Hankel function replaced by the simple power-law, the definite integration can be 
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done for each term in this series expansion. And it turns out that we can finally re-sum this 
series and get a closed- form result. 
To do this we rewrite 



HfXz) = J ilt (z) + iY^(z) = [1 + C0th(7T//)] J ifx (z) - 



sinh(7r/i) 



(4.5) 



The Bessel J, ltl function, together with some Xi-dependent coefficients, can be expanded as 



x- 1 1/2 J^{x 1 )e~^ = £ 



c+(xx) 



-i) n x 
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n\^T(n + 1 + 2i/i) 



(4.6) 



As noted, c n (xi) is proportional to a simple power of xi, thus can be integrate out in the x\ 
integration. The result is ^ 2 

7 + -- 



dxi X(x 1 )ct(x 1 ) 



4(-l) n e 



7r(2fi — i — 2m) 2 
The function j£ can be re-summed with respect to n, 
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where 



V^ (1) (^) = 



d ( dT{z)/dz 
d~z 1 T(z) 



(4.9) 



Similarly, the Bessel J_j„ function, together with its Xi-dependent coefficients, can be ex- 



1 Direct calculation, after an asymptotic expansion of the generalized hyper-geometric functions, shows 

4 e (At-m)7T 



7r(2/i — i — 2in) 



where /« ' are independent of x\. Note that, in (4.3), the hyper-geometric function contains an oscillation 
factor e~ 2lXl at x\ — > oo. This behavior makes either convergent or contain the similar oscillatory 
behavior, after we evaluate the result of the indefinite integral at x\ — > oo. The above expansion contains 
finite number of terms for each n and does not have this kind of oscillatory behavior. So it must be finite. 
This shows that all /„" must vanish, leaving only the first finite term. Explicitly we have only checked two 
of them, fn +n ' + ^ anc j 1 / 2 +«>+) ^ m( jeed van ish. Similar conclusion holds for J~ later. 

2 It is useful to deform the factor x 1 in (4.4 1 to x 1 1 ' 2+a by a small number a. This avoids some 
spurious divergences coming from some hypergeometric and Gamma functions in the intermediate steps. 
After the integration is done, we can take a — > limit and get a regular answer. 
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panded as 



2»-*f(_i)»Xi 
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The Xi integration can be performed term by term in n as 

4(-lV t 
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We re-sum over n and get 

oo 
n=0 

Combining the results we get 
C 2 = 
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A plot of C2 is shown in Fig. [3j again together with the numerical check. 
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(4.11) 



(4.12) 



(4.13) 



In the large M/H limit, the second line of (4.13) is exponentially suppressed, and the 



first line is power-law suppressed thus dominates. As a result, 

hm Co = 7: ■ 

^00 AM 2 



(4.14) 



Thus there is no Boltzmann factor in C 2 - in the large M limit, C2 is much greater than C\, 
and is the dominate contribution in this limit. 



5 Analytical continuation to imaginary fi 

In this section we extend our results to the case of lighter isocurvaton, < M < 3H/2, 
by analytically continuing /i — > —iv. As we have seen in [lj|2], both C\ and C 2 should have 
spurious divergences at u = 1/2, but those divergences cancel each other. These properties 



precisely show up after the analytical continuation of (3.7) and (4.13), because cosh (x) in 
(3.7) has a pole at x — —in/2 and ip^(x) in (4.13) has a pole at x = 0. The final result 



also matches well with the numerical results in Fig. 6 of Ref . |2] . These are shown in Fig. [4 
here. This confirmation provides a non-trivial check that we indeed have the exact results. 
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Figure 4: The analytical continuation of Ci(/i) and C2(a0 to imaginary /i (corresponds to 
real < v < 3/2). The solid line is C = C\ + C2. The dots are numerical results from [2]. 

6 An IR expansion for large /i 



In the previous sections, we have obtained the power spectrum correction from isocurvaton 
with arbitrary mass in quasi-single field inflation. However, if one is interested in more gen- 
eral scenarios than the present model, the exact method in the above sections could become 
very difficult or impossible to apply. Thus in this section, we present an approximation 
method - the IR expansion method. This method is applicable only in the large \x limit, 
and when the results are power-law (instead of exponential) suppressed, as we will explain. 
But on the other hand it involves simpler special functions and may be more friendly to be 
applied in more complicated cases. 

In the small argument limit, the Hankel function can be expanded as 



1 + COth(7T/i) 

r(i + i A t) 



7T 



-l /A 



(6.1) 



For small x, the expansion parameter for the Hankel function is ~ 4 ( n _j! 1+M ) 
So this leading order approximation is only good for < \x 



for terms of 
< y/Ji. We 



order 0(x 2n±i » 

emphasize that physically this is more than a good approximation in the IR limit. It includes 
the behaviour of the massive mode function at the super-Hubble limit kr — > 0, as well as an 
important portion of sub-Hubble scales —^JM/H < kr < 0. The behavior for larger x, 



H$(x) 



-1) 



3/4 




TlX 



X > jl 



(6.2) 



is not well represented at all by (6.1). The approximation method here is to replace the 



correct deep UV behavior in (6.2) with the UV behavior in (6.1 ), and run the integrals from 
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to oo with the approximation (6.1). 



To use this method, we need to satisfy the following two requirements for the integrals of 
interest. After a Wick rotation x = ±iy (which is explained below), we refer to y < JJL as 
the IR region, and y > JJi as the UV region. First, the UV behavior of (6.1 ) only contributes 



to the final result a term that is exponentially suppressed by large \i\ second, the correct UV 
behavior of the Hankel function also gives an exponentially small contribution. Then, if this 
approximation method gives a result that is power-law suppressed (instead of exponentially 
suppressed) by large //, this result is the correct leading order approximation. 

We demonstrate how this method works for the two integrals that we computed previ- 



ously. We first estimate the integral (3.1), which is schematically 



/•OO 

/ dyy- l ' 2 H^{iy)e 
Jo 



(6-3) 



Here it is important to note that we have made a Wick rotation x = iy. This rotation turns 
the oscillatory factor e tx into the exponential decay factor. Otherwise the UV contributions 
are not under as good control in the below estimates. 

First, we look at the UV contribution from the IR expansion method. In the large fi 



limit, (6.1) is schematically 
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Note that 



dy y 



-l/2±i Me -l/ _ e - 



TT/i/2 



(6.4) 



(6.5) 



So for yfji < y < oo, the integral (|6.3|) with the replacement (|6.4) contributes no more than 

(6.6) 



g — ITJt/2 _|_ g — 7T/V2 



~ e 
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Second, we estimate the actual UV contribution using (6.2). The approximation (6.2) is 
good for y > fi 2 , but can be used for y > ^JJl for our purpose of estimating the order of 
magnitude. The exponential factor e~ 2y in the integrand gives the suppression factor e -2 ^. 
So the actual UV contribution is of order 



~ e 



-VP 



(6.7) 



Finally, using the IR expansion method, plugging (6.1) in (6.3) and performing the resulting 
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integration precisely, we get 



8(1 + sin(2/iln2))-e" 2 ^ (6. 



in the large /i limit. 

To summarize the situation so far, we see that, although the UV contribution from both 
the IR expansion method and the actual integration are exponentially suppressed by large /x, 
the IR expansion method itself also gives an exponentially suppressed result. So we cannot 
get the exact leading order result using this method. However, we do know from the above 
estimates that C\ is exponentially suppressed. Thus it is negligible comparing with C 2 > as we 
have shown in previous sections, and shall show again below. 



We now look at the more important integral (4.1), which is schematically 



/*oo /»oo 

~e-*M d yi y^ /2 H^(-i yi )e^ dy 2 y~ 1/2 H%(-iy 2 ) e"" . (6.9) 

J0 Jyi 

Again here we have Wick-rotated such that X\ = —iyi, x 2 = —iy 2 - The direction of Wick 
rotation here is different from that we have used for C\. This is because we need to keep the 
exponential factor to be e~ yi and e~ V2 , specified by the ie-prescription of the vacuum choice 
in the interacting field theory. 

First, we estimate the UV contribution from the IR expansion method. After expanding 



the Hankel function using (6.1), for yi > // the inner integration gives, schematically, 

~ (vT + Vx") Vi 1/2 e- yi ■ (6.10) 



The above estimation can be extended to the region y > JJi for our purpose. Combined 
with the integrand in the outer layer, the leading UV contribution is 



~ e 



/•oo 

/ dy x y^e^e- 2yi ~ e" 2 ^ . (6.11) 
Second, we estimate the actual UV contribution. Similarly as we did previously, we use 



(6.2) and extend it to y > JJi for our order of magnitude purpose. Due to the exponentially 



suppressed factor, the UV contribution is of order 

~ e" 2 ^ . (6.12) 

To summarize, we see again that the UV contributions from both the IR expansion 
method and the actual integration are exponentially suppressed. In the following we will 



precisely compute the entire integration (4.1 ) using the IR expansion method (6.1 ), and show 



that it gives a power-law suppressed, hence the correct leading order result. 



We replace both Hankel functions in (4.1) with (6.1). The inner layer can be directly 
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computed and gives 



(-l) 1/ \2i)- 1 ^ 
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nT(-ifx) 
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However, combined with the integrand in the outer layer, the indefinite integral cannot be 
performed by Mathematica 8. We again use the trick of resummation - expanding the outer 
integrand in terms of power series, performing the simpler integration with the power-law 
terms, and re-summing the series. The power-law expansion of the integrand in the outer 
layer is 



n=0 



(-l) 1/4 (2z)-^ 



1 + COth(7T/i) (-l) n -1/2+^+ 



+ (-l)- 1/4 (2z)^ 
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T(l + ifi) n 

With the power-law, the final indefinite integral is doable using the following formulae 



(6.14) 



/x a 
dx x~ 1+a E u (ax) = [E u (ax) — Ei_ a (ax)] , 

a + v — 1 



(6.15) 



where 



EJz) 



oo e - zt 



dt = z v - l T(l - u,z) 



(6.16) 



is the exponential integral function. The resummation using the coefficients in (6.14) can 
be done by Mathematica. The final result is long and we do not list it here. We are only 
interested in the large \x limit, which is 
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(6.17) 



the same as (4.14). 



7 Truncating the interaction history 

As we have seen in the case of < M/H < 3/2 [Tj,[2] , the massive mode contributes to the 
curvature perturbation at different stages in its evolutionary history, sensitively depending 
on its mass. The lighter the isocurvaton is, the longer the contribution lasts after Hubble 
crossing. This is crucial in determining the shape of non-Gaussianities. 

To investigate where the dominant contribution comes from for the heavier mode M/H > 



3/2, we fix the upper limit of the integrals (3.1 ) and (4.1 ) to be oo and vary the lower limit 
x c 



3 Note that the IR cutoff x c is different from the cutoff we considered in Sec. |6] In Sec. |6j the cutoff (as a 
separation between UV and IR regimes) is taken after performing a Wick rotation. Thus it is completely a 
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There are at least two purposes for doing this: First, it is interesting to know the con- 
version from isocurvature to curvature perturbations is dominated by which period of time. 
Second, one could imagine models in which the turning of trajectory ends at some time 
before reheating. We hope this truncation of interaction can illustrate part of this effect]^] 

We numerically compute the integrals as a function of x c for different \x and plot them 
in Fig. |5j The numerical result shows that C2 decays quickly when we push x c into the far 
UV regime. On the other hand, it is interesting to see that the maximum UV contribution 
of C\ is actually also ~ 1/fi 2 for large //, but the effect of C\ decays a few e-folds after the 
Hubble crossing. As a result, for large //, the net contribution of C\ and C 2 does not grow 
monotonically in time (although the effects are small). This (if other effects do not wash 
it away) is a counter example of a no-go theorem stating that the curvature perturbation 
cannot be suppressed by the isocurvature perturbation (see PP/l for a discussion of super- 



Hubble perturbations and 18 for both sub- and super- Hubble). The no-go theorem is 



violated here because the isocurvature direction does not satisfy the slow-roll conditions, 



which is an assumption in 18 



In the case of a finite turn trajectory, Fig. [5] can also be reinterpreted as the contribution 
of the massive isocurvaton to the power spectrum as a function of comoving wave number 
k, with the replacement x c — > k/(a c H). In this way the pattern in Fig. [5] is a contribution 
to the power spectrum as a function of k due to the finite turn. 

Overall, for the constant turn case we are interested in in the previous sections, during 
the whole evolution, the massive isocurvaton contributes most significantly during the time 
interval from of order In fi e-folds before the Hubble crossing, to a few e-folds after it. Several 
e-folds after the Hubble crossing, the amplitude of the massive mode decays away. However it 
is a special feature of the massive mode with /i>l that starting from In [i e-folds before the 
Hubble crossing, the isocurvaton contribution already takes the similar order of magnitude 
as its eventual value. This feature can be observed in Fig. [6j The Ci, which is negligible after 
the Hubble-crossing, now becomes the dominate contribution before the Hubble-crossing, for 
\i 3> 1. Using results in Sec. |3j this term can be studied both numerically and analytically 
even with the cutoff x c . It has a maximal value « 0.05//i 2 , at x c ~ 0.5/x. 

Because the dominant isocurvature contribution comes from within a few e-folds, the 
constant-turn assumption made so far is not essential for the main results of this paper. 



mathematical trick, which does not change the hnal result of the integration but does change the intermediate 
behavior of the integration. While here, we impose a cutoff in real time instead of imaginary time. Thus 
the cutoff x c here has the physical meaning as the termination time of interactions. To keep the x c as the 
real-time cutoff, we can only Wick- rotate after shifting x — > x + x c , as we describe in Appendix [A| 

4 When the turning of trajectory ends, the physical process could be complicated: (1) the transfer from 
isocurvature to curvature terminates; (2) there may be classical oscillation because of the inertia of the 
inflaton; (3) the potential which stops the turning of trajectory may have extra features; (4) part of the 
vacuum state of the isocurvature direction is projected into a non-vacuum state of the inflaton direction; 
. . . . Here by putting a cutoff x c , we only included (1). Nevertheless this contributes part of the observable 
effects. Moreover, (2) and (3) could be eliminated by fine-tuning parameters in the potential, in which cases 
those effects may be separated. 
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Figure 5: The coefficients C, C\ and C 2 as a function of IR cutoff x c . The left, middle, right 
panels are with \i = 1, /z = 5 and /i = 10 respectively. 
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Figure 6: When an IR cutoff x c is considered, C\ has a maximum when considered as a 
function of x c (see the plot of C\ in Fig. |5|. The left panel plots the peak position x c which 
maximizes Ci, for different /i. The right panel plots the corresponding maximum values of 
C\ for those fi values. For large fi, one can read from the plot that the peak position is at 
x c 0.5/i, where C\ ps 0.05//1 2 . 
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As long as the variation of various parameters satisfy the conditions of the slow-variation 
approximation, for example, 8/(H6) <C 1 and M/(HM) <C 1, the same results apply with 
all the parameters take their instantaneous values. 



Notes added: When this paper was in completion, we became aware of a related work 19 
We thank the authors for coordinating on the submission of both works. 
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A Numerical computation of the integrals 

Although we have got full analytic results in Sec. [3]&[4| it is useful to check the calculation 
numerically. Note that the integrands in both integrals are oscillating at x\ — > oo. The 
oscillation makes the integration difficult (though still possible) to perform. There are a 
few techniques available to deal with such oscillations. For example, one can insert an e~ exi 
term, keeping e to be a small number. Or one can integrate the expression by parts to obtain 
better convergence. Here we use a third method, as used in (21: We Wick rotate X\ and x 2 
into imaginary values, such that e ±lXl , e ±iX2 — > e~ Xl , e~ X2 . Here the directions of rotation 
are chosen to ensure the convergence of the integrals, corresponding to the ie prescription in 
interacting field theories. 

For example, in C 2 , we first shift x\ — > X\ + x c , X2 — > x 2 + X\ and then rotate Xi — > —iyi, 
x 2 —> —iy 2 - We have 

/»oo /•oo /»oo /*oo 

/ dx\ I dx 2 f(xx,x 2 ) = / dxi / dx 2 f(xx + x c ,x 2 + xi) 

Jx c J xi JO JO 

roo /*oo 

= -/ dyi / dy 2 f(-iy 1 + x c ,-iy 2 -iy 1 ) . (A.l) 
Jo Jo 

Here / denotes the integrand in C 2 , and x c is an IR cutoff for the integral, added for two 
purposes: (1) In Fig. [5j we have plotted C 2 as a function of x c . Note that in this case, the 
Wick-rotation after the shifts does not change the fact that x c is the real time. (2) Even 
in the case without an IR cutoff, we actually also put a small x c in the above integral to 
prevent evaluating into 0/0 in the integrand. In the latter case we have checked that the x c 
we put (x c = 10 -10 ) is small enough such that the result is stable against order-of-magnitude 
variation of x c . 
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